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Abstract. The status of boson stars as black hole mimickers is presented among other 
mimickers. We focus on the analysis of the emission spectrum of a simple accretion disk model. 
We describe the free parameters that allow a boson star to become a black hole mimicker and 
present an example of a particular astrophysical case. 



1. Introduction 

One of the most important problems in relativistic astrophysics is that of the nature of black 
hole candidates (BHCs). It is commonly assumed that black hole solutions are correspond to 
BHCs. Nevertheless there are various other types of solutions of Einstein's equations sourced by 
different matter models that provide similar behavior to high energy observations in the vicinity 
of BHCs. These solutions are known as black hole mimickers. Among those solutions there 
are wormholes [U [2] , gavastars [21 El HI E] , brane world solutions [B] , naked singularities [7] and 
boson stars [SJ El UHl HI]- In fact, one important issue is to determine whether or not BHCs 
have an event horizon [K]. 

In order to estimate observable differences between black holes and mimickers is focused on 
the one hand to the search of gravitational wave fingerprints of the BHCs and their perturbative 
properties. In the first place perturbations to Schwarzschild black holes are very well studied (for 
a recent review see e.g. [I3]). On the other hand, black hole mimickers' perturbations have not 
been studied as much, however some results are in turn, for instance, the stability of the whole 
set of wormhole solutions usually thought of as mimickers is not clear, instead the instability 
of wormhole solutions has been shown for basic wormhole solutions \14:\ \T5\ [T6l I17j . destroying 
previous hope on the possibility that these are stable as shown in [18] for particular types of 
perturbations. Concerning gravastars, stability has been explored and has been found that there 
are stability regions [3l [El [20]. About boson star solutions the stability of the solutions has 
actually been exhaustively studied, for instance using perturbative methods [211 [22], catastrophe 
theory [23] and full non-linear numerical relativity, both in spherical symmetry [221 |25| and 
full 3D [26], and the stable branches of solutions are well known, and it has also been studied 
the fundamental quasinormal modes and tails of boson stars [27j . 

The study of boson stars has been pushed forward up to the binary boson star collision 
[28j and the non-linear evolution of perturbed boson stars [2^, in both cases considering the 
system as a source of gravitational waves. The study of gravitational wave signatures has also 
been studied in order to tell between a black hole and a gravastar Instead, in the case 



of wormholes e.g. simple solutions supported by a phantom scalar field would not stand the 
perturbative analysis and there is no hope for a study of a binary system, because as shown 
in [16] the life time of these solutions is rather short and they should either collapse and form 
black holes or explode before they could merge. 

On the other hand, a yet more common strategy to find observable differences among BHCs is 
associated to the study of the electromagnetic -observable- emission spectrum of geometrically 
thin optically thick disk models, since they are based on the study of time-like geodesies on 
a fixed background space-time and assume very simple conditions on the accreting material 

[a El El El [ini El. 

2. Boson stars 

Boson stars (BSs) are spherically symmetric stationary solutions to Einstein's field equations 
minimally coupled to a complex scalar field sourced by the stress energy tensor T^i, = 
l[df,^*d,(l) + df,(f>d,^*] - lgf,u[<l>*'''<P,a + V{\(l)\% with V{\<P\) = m2|0|2 + A|./-|4 where m is 
associated to the mass of the boson and A to the self-interaction of the boson system and 
the scalar field obeys the Klein-Gordon equation. The condition to find BSs solutions is that 
the scalar field has a harmonic dependence on time (j){r,t) = (j)Q{r)e~^^^ , where r is the radial 
coordinate, t the coordinate time and a; is a characteristic frequency that determines the central 
density of the solution [10]. This condition implies that the stress-energy tensor becomes time- 
independent, which in turn implies that the geometry of the space-time is also time- independent. 
Using the rescaling c^o = V^(t>o, f = mr, i = cot, a = and A = g^^i , one constructs the 
boson star solutions like those shown in Fig la. for various values of A. Each point of each curve 
represents a BS solution. Those to the left of the maxima are stable whereas those to the right 
are unstable [30] . 

3. The accretion disk model 

The disk model around a given BHC is defined as follows. It corresponds to a geometrically 
thin, optically thick, steady accretion disk. The power per unit area generated by such a disk 
rotating around a central object is given by [31] and used also in the study of disks around 

gravastars [5] and naked singularities [7] reads D{r) = ■Spai~w) {e-ql)'^ IrS^ ~ ^^)^^'^' 
where M is the accretion mass rate, E, L, are the energy per unit mass, angular momentum 
per unit mass and angular velocity of test particles the disk is made of and rj is the inner edge 
of the disk. For black holes this radius is assumed to be at the ISCO (r = 6M) of the hole. 
For BSs we choose rj = based on two considerations: i) BSs allow circular orbits in the whole 
spatial domain and thus there is no kinematical restriction to choose = and ii) it has been 
shown that the luminosity of the disk for BSs considering this inner radius never reaches the 
value of Eddington luminosity in the whole spatial domain even for high accretion rates and 
thus there are no radiation pressure effects imposing a restriction on the inner edge location 
[9]. Now, assuming it is possible to define a local temperature the Steffan-Boltzmann law is 
used so that D{r) = aT'^, where a = 5.67 x 10~^ erg cm~'^ K~'^ is the Steffan-Boltzmann 
constant. Now, considering the disk emits as a black body, we use the dependence of T on 
the radial coordinate and therefore the luminosity L{i') of the disk for each frequency i' and 
the fiux F{iy) can be calculated using the expression for the black body spectral distribution: 
L(z>) = 4TT(fF{i') = cos{'d)u^ J^^ ^hJ/kT_i 1 where d is the distance to the source, and 
rf indicate the location of the inner and outer edges of the disk, h = 6.6256 x 10"^'^ erg s is 
the Planck constant, k = 1.3805 x 10~^^ erg is the Boltzmann constant and is the disk 
inclination. The algorithm to construct the emission spectrum for such a model of accretion 
disk around a BS and a BH is as follows: 1) Define the space-time functions a and a by choosing 
one of the equilibrium configurations in Fig. la and calculate M; 2) Define the metric of 



the equivalent BH through ubh = y^l — 2M/r and gbh = 1/a/1 — 2M/r; 3) Calculate the 
angular velocity, angular momentum and energy of a test particle for both space-times 0,bs,bh, 
Lbs,bh, Ebs,bh', 4) Use such quantities to calculate the power emitted in both cases DBs{r) 
and DBH{r)', 5) Calculate the temperature of the disk in both cases TBs{f) = {DBsiT) / ^Y^^ 
and TBH{f) = {DBH{f) / crY^^', 6) Use such temperature to integrate the luminosity Lbs{'^) and 
Lbh{^) for several values of v. 



4. Mimicker configurations 

The free parameters of BSs that one can play with in order to produce a spectrum that mimics 
that due to the presence of a black hole are m and A. As an example in Fig. lb the power 
spectrum of a black hole is compared with two spectra due to two different BS configurations of 
the same mass but with different self- interaction A, one of which is a mimicker. It can be noticed 
how the spectrum due to one of the BSs lies on top of that due to the black hole. This is the 
effect that a black hole mimicker produces. For each value of A it is possible to find mimickers 
by changing m. For instance, in Fig. la stars indicate the configurations corresponding to black 
hole mimickers for various values of A. 




Figure 1. In the left panel it is shown a set of BSs configurations; each point of each curve 
corresponds to a BS configuration. Those configurations to the left of the maxima are stable, 
those between the maxima and the The mimicker configurations are indicated with stars for 
values of A = 0, 10, 20, 30, 40, 50, 60, 70. It is expected that for every value of A there is a 
mimicker configuration that lies near the ones indicated for the values shown here. We show the 
mass of the inverted triangles are unstable with negative binding energy and collapse into black 
holes, whereas those to the right of the inverted triangles are unstable configurations that instead 
of collapsing explode |30j . In the right panel we show an example of spectrum generated by a disk 
around BS and around a Schwarzschild black hole. The mimicker configuration used corresponds 
to a BS configuration with M = 0.9898Mp;/m, A = 50 and boson mass m = 5.25 x 10~^^GeV, 
and works as a mimicker of a black hole with mass M = IOMq with an accretion rate parameter 
2 X 10~^^MQ/?/r. The spectra due to the black hole and to the mimicker are on top of each 
other, and for reference the spectrum due to a non-mimicker configuration is also shown. 



5. Discussion 

A summary of boson stars as black hole mimickers has been presented. It has been found that 
it is possible to find a boson star configuration that may supplant a black hole based on the 
emission spectrum of an accretion disk model. Nevertheless there are items that have to be ad- 
dressed if boson stars are to still be considered black hole mimickers. The most important one 



is the question of whether or not boson stars are capable of trapping baryons, which eventually 
would form a mechanical surface that would behave like that of a baryonic compact object as as- 
sumed in studies of mimickers [12] . The fact that boson star space-times are regular everywhere 
indicates that particles (e.g. baryons) following geodesies would truly be accreted and retained 
on top of the boson star space-time; this would imply that there is nothing like accretion of 
matter, and instead, that the region of the boson star is a place where particles fall into and are 
expelled out, because it has been shown [TU] that particles following radial trajectories would 
cross the origin and get out, and when test particles have a non-zero impact parameter there 
is an infinite potential barrier that forbids test particles to remain near the center of the boson 
star. However this is not a definitive answer, instead it is necessary to consider the self-gravity 
of baryons and their contribution to the geometry. With that information at hand it will be 
possible to determine whether or not boson stars are still potential black hole mimickers and 
then proceed to other more stringent astrophysics related tests like gravitational lensing and 
gravitational radiation. 
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